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ABSTRACT 

Aims. Although the time of the Maunder minimum (1645-1715) is widely known as a period of extremely low solar activity, claims are still 
debated that solar activity during that period might still have been moderate, even higher than the current solar cycle # 24. We have revisited all 
the existing pieces of evidence and datasets, both direct and indirect, to assess the level of solar activity during the Maunder minimum. 
Methods. We discuss the East Asian naked-eye sunspot observations, the telescopic solar observations, the fraction of sunspot active days, the 
latitudinal extent of sunspot positions, auroral sightings at high latitudes, cosmogenic radionuclide data as well as solar eclipse observations for 
that period. We also consider peculiar features of the Sun (very strong hemispheric asymmetry of sunspot location, unusual differential rotation 
and the lack of the K-corona) that imply a special mode of solar activity during the Maunder minimum. 

Results. The level of solar activity during the Maunder minimum is reassessed on the basis of all available data sets. 

Conclusions. We conclude that solar activity was indeed at an exceptionally low level during the Maunder minimum. Although the exact level 
is still unclear, it was definitely below that during the Dalton minimum around 1800 and significantly below that of the current solar cycle # 24. 
Claims of a moderate-to-high level of solar activity during the Maunder minimum are rejected at a high confidence level. 

Key words. Sumactivity - Sun:dynamo - Sun:Maunder Minimum 


1. Introduction 

In addition to the dominant 11-year Schwabe cycle, solar ac¬ 
tivity varies on the centennial time scale ( Hathaway! Il2010j) . 
It is a common present-day paradigm that the Maunder min¬ 
imum (MM), occurring during the interval 1645-1715 ( Eddy , 
1976h . was a period of greatly suppressed solar activity called 
a Grand minimum. Grand minima are usually considered as 
periods of greatly suppressed solar activity corresponding to a 
special state of the solar dynamo dCharbonneaui [201 Oh . Of spe¬ 
cial interest is the so-called core MM ( 1645-1700) when cycli c 
sunspot activity was hardly visible ( Vaquero & Trigol 120151) . 
While such Grand minima are known, from the indirect evi¬ 
dence provided by the cosmogenic isotopes 14 C and 10 Be data 


for the Holocene, to occur sporadically, with the Sun spend 


ing on average 1 /§ of the time in such a state ( Usoskin et al 


20071) . the MM is the only Grand minimum covered by di¬ 


rect solar (and some relevant terrestrial) observations. It there¬ 
fore forms a benchmark for other Grand minima. We note 
that other periods of reduced activity during the last cen¬ 
turies, such as the Dalton minimum at the turn of the 18th 
and 19th centuries, the Gleissberg minimum around 1900, or 
the weak present solar cycle #24, are also known but they are 


typically not consider ed to be Grand minima ( Schiissler et al 


1997; Sokoloff, 2004). However, the exact level of solar ac¬ 


tivity in the 17th century remains somewhat uncertain (e.g . 
Vaquero & Vazque j. 120091: IVaquero et all 1201 ll: Idette et al 


20141) . leaving room for discussions and speculations. For ex- 
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ample, there have been several suggestions that sunspot activity 
was moderate or even high during the core MM (1645-1700), 


cle #24 (SchovjJ 

9§I Gleissberg et al] 1979; Culler] 1980; 

Nagovitsyn, 1997 

Ogurtsov et al.. 2003|: Nagovitsyn et al. 

2004 

; Volobuev, 

2004; Rek, 20 FI Zolotova & Ponvavin 

2015 

). Some of these were based on a mathematical synthe- 


sis using empirical ru les in a way similar to ISchovel (119551) 
and INagovitsvnl ( 1997 ). and therefore are not true reconstruc¬ 
tions. Some others used re -analysis of the direct data series 
jRekl 120131: [Zolotova & Ponvavin . 12015 ) and provide claimed 
assessments of the solar variability. While earlier suggestions 
have been convincingly rebutted bv lEddvl ( 1983 ), the most re¬ 
cent ones are still circulating. If such claims were true, then 
the MM would not be a Grand minimum. This would poten¬ 
tially cast doubts upon the existence of any Grand minimum, 
including those reconstructed from cosmogenic isotopes. 

There are indications that the under lying solar mag 


netic cycles still op erated during the MM ( Beer et all 1 19981: 


lUsoskin et a U 2£pi^but at the threshold level as proposed al¬ 
ready by lMaunder (1922 


It ought not to be overlooked that, prolonged as this in¬ 
activity of the Sun certainly was, yet few stray spots 
noted during “the seventy years’ death”- 1660, 1671, 
1684, 1695, 1707, 1718 [we are, however, less certain 
about the exact timings of these activity maxima]- cor¬ 
respond, as nearly as we can expect, to the theoretical 
dates of maximum. ... If I may repeat the simile which 
I used in my paper for Knowledge in 1894, “just as in 
a deeply inundated country, the loftiest objects will still 
raise their heads above the flood, and a spire here, a hill, 
a tower, a tree there, enable one to trace out the con¬ 
figuration of the submerged champaign,” to the above 
mentioned years seem be marked out as the crests of a 
sunken spot-curve. 


The nature of the MM is of much more than purely 
academic interest. A recent analysis of cosmogenic isotope 
data revealed a 10% chance that Maunder minimum condi¬ 
tions would return within 5 0 years of now (Lockwood . 2010t 


Solanki & Krivoval 1201 lk iBarnard et all 1201 ll) . It therefore be¬ 
comes of practical importance to accurately describe and un¬ 
derstand the MM, since a future Grand minimum is expected 
to have significant implications for space climate and space 
weather. 


Here we present a compilation of observational and his¬ 
torical facts and evidence showing that the MM was indeed a 
Grand minimum of solar activity and that the level of solar ac¬ 
tivity was very low, much lower than that during the Dalton 
minimum as well as the present cycle # 24 . In Section[2]we re¬ 
visit sunspot observations during the MM. In Section [3] we an¬ 
alyze indirect proxy records of solar activity, specifically auro¬ 
rae borealis and cosmogenic isotopes. In Section |4] we discuss 
consequences of the MM for solar dynamo and solar irradiance 
modelling. Conclusions are presented in Section 0 



Fig. 1 . Annual group sunspot num bers during and around t he 
Maunder minimum, according to Hoyt & Schatte n (1998) - 
GSN, IZolotova & Ponvavinl d2015l) - ZP15, and loose and 
strictly conservative models from Vaquero et al.1 ( 2015al) (see 
Sect. 12. IK as denoted in the legend. 


2. Sunspot observation in the 17th century 

Figure |T] shows different estimates of sunspot activity, quan¬ 
tified in terms of the annual group sunspot number (GSN) 
/?G,_around the MM. The conventional GSN dHovt & Schatteri 


19981 called HS98 henceforth), with the recent corrections, re¬ 


lated to newly uncover ed data or connections of earlier errors , 
applied (see details in Vaguero et all 1201 lb I Vaquero & Trigol 


l20l4 ILockwood et al. . 2014bl) . is shown as the black curve. 


This series, however, contains a large number of generic no¬ 
spot statements (i.e. statements that no spots were seen on the 


tion (Kovaltsov et al., 2004; 

Vaquero, 2007; Clette et al., 2014; 

Zolotova & Ponvavin, 2015 

; Vaquero et al., 20154 see also 


Sect. 12.31) . Figure [I] shows also two recent estimates of the an¬ 
nual GSN by IVaquero et al.1 d2015al) . who treated generic no¬ 
sunspot records in the HS98 catalog in a conservative way. The 
sunspot numbers were estimated using the acti ve-vs -inactive 
day st atistics (see Sect. 12.11 full details in IVaquero et al 


2015al) . All these results lie close to each other and imply 


very low sunspot activity during the MM. On the contrary, 


Zolotova & Ponyavin (120151) . called henceforth ZP15, argue 


for higher sunspot activity in the MM (the red dotted curve 
in Figure |T] is taken from Figure 13 of ZP15), with the sunspot 
cycles being not smaller than a GSN of 30 and even reaching 
90-100 during the core MM. 

For subsequent analysis we consider two scenarios of so¬ 
lar activity, reflecting the opposing views on the level of 
solar activity around the MM before 1749: (1) L-scenario 
of low activity during the MM, as based on the conven¬ 
tional GSN ( iHovt & Schattenll 1998) with the recent correc¬ 
tions implemented (see ILockwood et all l2014bi for details) 
- see black curve in Fig. \T\ (2) H-scenario of high activ¬ 
ity during the MM, based on GSN as proposed by ZP15 
(red dotted curve in Fig. []]). This scenario qualitatively repre¬ 
sents also other sugge s tions o f high activ i ty (e.g ., INagovitsvnl 
19971: lOgurtsov et all 120031: IVolobuevi 12004 ). After 1749, 


both scenarios are extended by the International sunspot num¬ 
ber (http://sidc.oma.be/silso/datafiles). We use annual values 
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throughout the paper unless another time resolution is explic¬ 
itly mentioned. 


2.1. Fraction of active days 

High solar cycles imply that ~ 100% of days are active, 
and sunspots are seen on the Sun almost every day dur¬ 
ing such cycles, except for a few years around cycl e min¬ 
ima dKovaltsov et all l2004 IVaquero et all 1201 2l 120141) . If the 


sunspot activity was high during the MM, as proposed by the 
H-scenario, the Sun must have been displaying sunspots almost 
every day. However, this clearly contradicts with the data, since 
the reported sunspot days, also those reported by active ob¬ 
servers, cover only a small fraction of the year even around the 
proposed cycle maxima (see Fig. 2 in IVaquero et all Eoi5ah . 
Thus, either one has to assume a severe selection bias for ob¬ 
servers reporting only a few sunspot days per year while spots 
were present all the time, or to accept that indeed spots were 
rare. 

During periods of weak solar activity, the fraction of 
spotless days is a very sensitive indicator of the ac 


tivity level (Harvey & White! 1 19991: iKovaltsov et all 12004 


IVaquero & Trig o. 2014), much more precise than the sunspot 


counts. However, this quantity tends to zero (almost all days 
are active) when the average sunspot num ber exceeds 20 
( Vaquero et al.U2015al) . IVaquero et al l ( 2015a ) considered sev¬ 
eral statistically conservative models to assess the sunspot 
number during the MM from the active day fraction. The 
’’loose” model ignores all generic no-spot statements and ac¬ 
cepts only explicit no-spot records with exact date and ex¬ 
plicit statements of no spots on the Sun, while the ’’strict” 
model considers only such explicit statements as in the ’’loose” 
model but made by at least two independent observers for the 
spotless days. In this way, the possibility of omitting spots is 
greatly reduced since the two observers would have to omit the 
same spot independently. The strict model can be considered 
as the most generous upper bound to sunspot activity during 
the MM. However, it most likely exaggerates the activity by 
over-suppressing records reporting no spots on the Sun. These 
models are shown in Fig. |T] One can see that these estimates 
yield sunspot numbers that do not exceed 5(15) for the ’’loose” 
(’’strict”) model during the MM. 


2.2. Occidental telescopic sunspot observations: 
Historical perspective 

The use of the telescope for astronomical observations be¬ 
came widespread quickly after 1609. We know that there were 
telescopes with sufficient quality and size to see even small 
spots, in the second half of the 17th century. It is also known 
that astronomers of that era used other devices in their rou¬ 
tine observations such as mural quadrants or meridian lines 
jHeilbornlll999l) . However, as proposed by ZP15, the quality 
of the sunspot data for that period might be compromised by 
non-scientific biases. 


2.2.1. Dominant world view 

Recently, ZP15 suggested that scientists of the 17th century 
might be influenced by the “dominant worldview of the sev¬ 
enteenth century that spots (Sun’s planets) are shadows from 
a transit of unknown celestial bodies”, and that “an object on 
the solar surface with an irregular shape or consisting of a set 
of small spots could have been omitted in a textual report be¬ 
cause it was impossible to recognize that this object is a celes¬ 
tial body”. This would suggest that professional astronomers of 
the 17th century, even if technically capable of observing spots, 
might distort the actual records for politically/religiously mo¬ 
tivated, non-scientific reasons. This was the key argument for 
ZP15 to propose the high solar activity during the MM. Below 
we discuss that, on the contrary, scientists of the 17th century 
were reporting sunspots quite objectively. 

Sunspots: Planets or solar features? 

There was a controversy in the first decades of the 17th century 
about location of sunspots: either on the Sun (like clouds), or 
orbiting at a distance (as a planet). However, already Scheiner 
and Hevelius plotted non-circular plots and showed the per¬ 
spective foreshorteni ng of spots near the limb. In his Accuratior 
Disquisito, Christoph Scheiner (1612) wrote pseudonymously 
as ‘Appelles waiting behind the picture’ and detailed the ap¬ 
pearance o f spots as of irregular shape and variable, and finally 
concluded (Gal ileo & Schemed 120 id) : 


They are not to be admitted among the number of stars, 
because they are of an irregular shape, because they 
change their shape, because they [... ] should already 
have returned several times, contrary to what has hap¬ 
pened, because spots frequently arise in the middle of 
the Sun that at ingress escaped sharp eyes, because 
sometimes some disappear before having finished their 
course. 

Even though Scheiner had believed until this point that 
sunspots were bodies or other entities just outside the Sun, he 
did note all their properties very objectively. Later, Scheiner 
( 1630 ) concluded in his comprehensive book on sunspots, 
“Maculae non sunt extra solem” (spots are not outside the 
Sun, p. 455ff.) and even “Nuclei Macularum sunt profundi” 
(the cores of sun spots are deep, p. 506). On the contrary, 
Smogulecz & Schonberger (1626) who were colleagues of 
Scheiner in Ingolstadt and Freiburg-im-Breisgau, respectively, 
called the spots “stellaae solares” (= solar stars) which was 
meant in the sense of moons. Some authors, especially anti- 
Copernican astronomers, suc h as Antonius Mari a Sch yrleus 
of Rheita (1604-1660) (see iGomez & Vaqueroi 12015 ) and 
Charles Malapert (1581-1630), followed the planetary model. 
On the other hand, Galileo had geometrically demonstrated 
(using the measured apparent velocities of crossing the solar 
disc) that spots are located on the solar surface. In fact, the 
changes in the trajectory of sunspots on the ’’solar surface” 
were an im portant element of discussion in the context of he¬ 
liocentrism dSmithl 1 19851: iHutchisoni Il990l: iTopped. Il999h . 

It was clear already in that time that sunspots are not plan¬ 
ets, for reasons of the form, color, shape of the spots near 
the limb and their occasional disappearance in the middle of 
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the disk. A nice example of the kind is given in a letter to 
William Gascoigne (1612-1644), that William Crabtree wrote 
on 7 Augu st 1640 (= Aug 17 greg.) dChapmani 120041) . as pub¬ 
lished bv lDerham & Crabtri9 (1171 ll) : 


I have often observed these Spots; yet from all my 
Observations cannot find one Argument to prove them 
other than fading Bodies. But that they are no Stars, but 
unconstant (in regard of their Generation) and irregu¬ 
lar Excrescences arising out of, or proceeding from the 
Sun’s Body, many things seem to me to make it more 
than probable. 


Although some astronomers still believed in the mid-17th 
century that sunspots were small planets orbiting the Sun, 
the common paradigm among the astronomers of that time 
was “that spots were current material features on the very 
surface of the Sun” (Brody . 120021. page 78). Therefore, ob¬ 
servers of sunspots during the MM, in particular professional 
astronomers, did not adhere to the ’’dominant worldview” of 
the planetary nature of sunspots and hence were not strongly 
influenced by it, contrary to the claim of ZP15. 

Galileo’s trial. 

We note that the problem in the trial of Galileo was not the 
Copernican system, but the claim that astronomical hypothe¬ 
ses can be validated or invalidated (an absurd presumption for 
many people of the early 17th c entu ry) and the potential claim 
of re-interpreting the Bible (iSchroderi 2002). The planetary 
system was considered as a mathematical tool to compute the 
motion of the planets as precisely as possible; it was not a sub¬ 
ject to be proven. This subtle difference was an important issue 
in the first half of the 17th century to comply with the require¬ 
ments of the catholic church. While an entire discussion of the 
various misconceptions about the Galileo trial is beyond the 
scope of this Paper, there are many indications that the nature 
and origin of celestial phenomena had been discussed by the 
scholars of the 17th century, rather than discounted by a stan¬ 
dard world view. We are not aware of any documental evidence 
that writing about sunspots was prohibited or generally disliked 
by the majority of observers in any document. 

Shape of sunspots. 

ZP15 presented a hand-picked selection of a few drawings to 
support their statement that “there was a tendency to draw 
sunspots as objects of a circularized form”, but there are plenty 
of other drawings from the same time showing sunspots of 
irregular shape and sunspot groups with complex structures. 
Here we show only a few of numerous examples. Figure [2] 
depicts a sunspot group observed in several observatories in 
Europe in August 1671. A dominant spot with a complex struc¬ 
ture having multiple umbrae within the same penumbra and a 
group of small spots in its surroundings can be appreciated. 
Another example (Figure [3]) shows a spot drawing by G.D. 
Cassini in 1671 (lOldenburgl Il67 lcb 0 , which illustrates the 
complexity and non-circularity as well as the foreshortening 


1 Henry Oldenburg was Secretary of the Royal Society and com¬ 
piled findings from letters of other scientists in the Philosophical 
Transactions in his own words. We therefore cite his name although it 
is not given for the actual article. 





Fig. 2. Drawing of a sunspot group observed in August 1671, 
as published in number 75 of the Philosophical Transactions, 
corresponding to August 14, 1671. 


of sunspots very clearly. Finally, Figure 0] displays a sunspot 
observation made by J. Cassini and Maraldi from Montpellier 
(Mar-29-1701). There is a small sunspot group (labelled as A) 
approximately in the middle of the solar disc which is zoomed 
in the bottom left corner, exhibiting a complex structure, and a 
legend that reads “Shape of the Spot observed with a large tele¬ 
scope”. These drawings are not limited to “circularized forms”, 
and such instances are numerous. 

It is important that observers who made drawings actually 
retained th e perspective foreshorten i ng of t he spots near th e 

~~ dl63f)h - iHevelinsI ill 647h 
1671c) as well as Cassini 


solar limb. Galilei 
G.D. Cassini in 1671 (lOldenbur; 
dl730l obser vation of 1 684), P. De La Hire! dl720i observation 
of 1703), and lDerhaml ( 1703 ) all drew slim, non-circular spots 
near the edge of the Sun. It was clear to them that those ob¬ 
jects cannot be spheres. They were not shadows either since 
that would require an additional light source similar to the Sun 
which is not observed. A note by G.D. Cassini of 1684 says 
dCassinil 1730): 


This penumbra is getting rounder when the spot ap¬ 
proaches the center, as it is always happening, this is 
an indication that this penumbra is flat, and that it looks 
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Fig. 3. Sunspot drawing of by G.D. Cassini in 1671 
dOldenburd. 1167 id) . 


narrow only because it is presenting itself in an oblique 
manner, as is the surface of the Sun towards the limb, 
on which it has to lie. 



Fig. 4. Sunspot observed by J. Cassini and Maraldi from 
Montpellier (Mar-29-1701). Reproduced from page 78 of the 
Histoire de L’Academie Roy ale des Sciences (Annee MDCCI). 


2.2.2. The very low activity during 1660-1671 

The period of 1660-1671 is indicative of very low activity in 
the HS98 database, but it is mostly based on generic statements 
of the absence of sunspots. For example, based on a report by 
G.D. Cassini, a sunspot observed in 1671 (lOldenburgi 1167 lb) 
was described in detail, and it was noted that 


“it is now about twenty years since, that Astronomers 
have not seen any considerable spots in the Sun, though 
before that time [... ] they have from time to time ob¬ 
served them. The Sun appeared all that while with an 
entire brightness.” 


While G.D. Cassini was an opponent of Copernicus and 
Newton (Habashi, 2007) and in fact discovered a number of 
satellites of Saturn, he did accept that sunspots are on the solar 
surface and did not alter their appearance to make them circu¬ 
lar. 


Thus, the idea suggested by ZP15 of the strong influence 
of theological or philosophical ideas about the perfection of 
the celestial bodies (especially the Sun) on professional as¬ 
tronomers in the late 17th century is not supported by our ac¬ 
tual knowledge of solar observations and scientific vision dur¬ 
ing that time. There is sufficient evidence that, beginning with 
the use of telescopes in astronomy, the existence and the na¬ 
ture of sunspots were thoroughly discussed, including various 
opinions, but based on the best observation technology of the 
time. We further conclude that sunspots were not omitted de¬ 
liberately from observing records for religious, philosophical 
or political reasons during the MM. The observational cover¬ 
age was just incomplete and partially vague. Moreover, many 
existing pieces of evidence imply that spots of different shapes 
were recorded, contrary to the claim of ZP15. 


The last sentence implies that the Sun was also void of any 
other dark features, even if they would not have been reported 
in terms of sunspots. There is also a footnote saying that in¬ 
deed some spots were witnessed in 1660 and 1661, so the 
20 years mentioned are exaggerated. The Journal also states 
(lOldenburgi. 11671 al) that “as far as we can learn, the last obser¬ 
vation in England of any Solar Spots, was made by our Noble 
Philosopher Mr. Boyl” on Apr 27 (= May 7 greg.) 1660 and 
May 25 (= Jun 4 greg.). He described a “very dark spot al¬ 
most of quadrangular form”. Additionally, one of the spots was 
described as oval, while another was oblong and curved. This 
statement contradicts the assumption by ZP15 that the majority 
of non-circular spots were omitted, especially in conjunction 
with the surprise with which the article was written that spots 
were seen at all. If there were a number of non-circular spots 
during this 10-year period (allegedly not reported), there would 
have been no reason to ‘celebrate’ yet another non-circular spot 
in 1671. 

As another example. ISporerl (1 889 ), p. 315, cited a note by 
Weigel from Jena in 1665 which can be translated as 
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Many diligent observers of the skies have wondered 
here that for such a long time no spots were notice¬ 
able on the Sun. And we need to admit here in Jena 
that, despite having tried in many ways, setting up large 
and small spotting scopes pointed to the Sun, we have 
not found such phenomena for a considerable amount 
of time. 


Since the notes on the absence of spots come from various 
countries and from catholic, protestant and Anglican people, 
we believe there was no wide-spread religious attitude to with¬ 
hold spots in order to save the purity of the Sun. 

The only positive sunspot report between 1660 and 1671 
in the HS98 database is the one by Kircher in 1667. This data 
point comes from a note (IFrickL 1681, p. 49) stating that 


the late Christoff Weickman, who was experienced in 
optics and made a number of excellent telescopes, 
watched the Sun at various times hoping to see the like 
[sunspots] on the Sun, but could never get a glimpse of 
them [... ] So Mr Weickman wrote to Father Kircher 
and uncovered him that he could not see such things on 
the Sun, does not know why this is or where the mis¬ 
take could be. Father Kircher answered from Rome on 
2 September 1667 that it happens very rarely that one 
could see the Sun as such; he had not seen it in such a 
manner more than once, namely Anno 1636. 


One can see that the date of the letter in 1667 was mistakenly 
considered as the observing date. Instead the report clearly in¬ 
dicates that no sunspots were seen at all by Weickman in the 
1660s. The sunspot observation by Kircher in 1667 is erro¬ 
neous and needs to be removed from the HS98 database. Then 
no indication of sunspots exists in the 1660s. We note that this 
false report was used by ZP15 to evaluate the sunspot cycle 
maximum around that date. 


2.3. Generic statements and gaps in the HS98 
database 

The database of HS98 forms a basis to many studies of sunspots 
records during the period under investigation. In particular, 
ZP15 based their arguments on this database without referring 
to the original records. However, it contains a number of not 
obvious features which can be easily misinterpreted if not con¬ 
sidered properly. Here we discuss such features which are di¬ 
rectly related to the evaluation of sunspot activity in the 17th 
century. 

In particular, many no-spot records were related to astro- 
metric observations of the Sun such as solar meridian a ltitude 
or the app arent solar diam eter ( Vaquero & Gallegol 12014 ). For 
example, iManfredil dl736l) listed more than 4200 solar merid¬ 
ian observations made by several scientists during 1655-1736 
using the gigantic camera obscura installed on the floor of the 
Basilica of San Petronio in Bologna. Such observations were 
not focused on sunspots and did not include any mentioning 
of spots. However, HS98 treated all these reports as observa¬ 
tions of the absence of sunspot groups which, of course, was 
incorrect. 


HS98 database contains gaps of the observing records of 
Marius and Riccioli, that occur exactly during days when other 
observers reported spots, which was interpreted by ZP15 as 
indications that they deliberately stopped reporting to hide 
sunspots: “It is noteworthy that when the Sun became active, 
Marius and Riccioli immediately stopped observations.”. We 
note that this interpretation is erroneous and based on igno¬ 
rance of the detail of the HS98 database as explained below. 

The original statement by Marius from Apr 16 (= Apr 26 
greg.) 1619, on which this series is based, is 

While I did not find as many spots in the disk of the 
Sun over the past one-and-a-half years, often not even a 
single spot, which was never seen in the year before, I 
noted in my observing diary: Mirum mihi videtur, adeo 
raras vel saepius nullas maculas in disco solis depre- 
hendi, quod ante hac nunque est observatum 


which is a Latin repetition of what he said before. Marius 
clearly states that the sunspot number was not exactly zero, but 
very low. HS98 have used this statement to approximate the 
activity by zeros in their database, more precisely by filling all 
dates of the 1.5-yr interval with zeros except the periods when 
other observers did see spots. The existence of the gaps is by 
no means based on the actual observing report by Marius, but 
is an effect of the way HS98 have interpreted the comment. 

The same reason holds for the gaps in the sunspots re¬ 
ported by Riccioli (1653), p. 96, whose data (zeros) in the HS98 
database are based on the statement that 


...in the year 1618 when a comet and tail shone, 
no spots were observed, said Argolus in Pandosion 
Sphaericum chapter 44. 


The original statement by Argolus (1644), p. 213, states: “Anno 
1618 tempore quo Trabs, et Cometa affulsit nulla visa est.” 
Apart from the fact that it was not Riccioli himself who ob¬ 
served, this again led to filling all days in 1618 with zeros (in 
the HS98 database) except the days when other observers saw 
spots. 

The method of filling the HS98 database for many months 
and even years with zeros is based on generic verbal reports on 
the absence of spots for long periods also in the cases of Picard, 
G.D. Cassini, Decha l es, Maraldi, Siverus a nd others (see, e.g, 
IVaquero et all 1201 ll: IVaquero et all l2015al) . HS98 must have 
filled those periods in the sense of probably very low activ¬ 
ity, but they are not meant to provide exact timings of obser¬ 
vations, as ZP15 interpreted them. The appearance of gaps in 
zero records when other observers reported spots is not an in¬ 
dication of withholding spots in observing reports but rather a 
simple technical way of avoiding conflicting data in the HS98 
database. ZP15 mistook the entries in the HS98 database for 
actual observing dates and interpreted them incorrectly. 

While assuming a large number of days without spots is 
a significant unde restimation of the sol ar activity on the one 
hand, as shown bv I Vaquero et al.1 ( 2015al) and has been pointed 
out by ZP15 as well, the assumption that observers deliberately 
stopped reporting is, on the other hand, not supported by any 
original text and remains an ungrounded speculation. 
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The observations by Hevelius of 1653-1684 as recovered 
bv lHovt & Schattenl (1995) should also be scrutinized with re¬ 
gard to a possible omission of spots. Citing the former refer¬ 
ence, ZP15 even claim that “Hevelius quite consciously did 
not record sunspots”, while the original statement was that 
“Hevelius occasionally missed sunspots but usually was a re¬ 
liable observer.” Actually, out of 24 groups that could have 
bee n detected by Hevelius given his observing days, he saw 
20 dHovt & Schattenl Il995l) . He never reported the absence of 
sunspots when others saw them. The four occasions are simply 
not accompanied by any statement about presence or absence 
of spots. This can be understood considering that the sunspot 
notes are just remarks on his solar elevation measurements 
iHeveliusl I 16791 part 3). Those, however, were made with a 
quadrans azimutalis which had no telescope, since Hevelius re¬ 
fused to switch to a telescope at some point, perhaps because he 
did not want to spoil his time series of measurements dHabashil 


2007). He therefore could not see sunspots at all with his de¬ 


vice and had to take an additional instrument to observe them, 
and it is probable that he did not do so on each day he mea¬ 
sured the solar elevation, hence left so many days neither with 
positive nor with negative information on sunspots. We have to 
treat those as non-observations. 


2.4. Methodological errors of ZP15 

The original work by ZP15 unfortunately contains a number of 
methodological errors which eventually led them to an extreme 
conclusion that sunspot activity during the MM was at a mod¬ 
erate to high level. In particular, ZP15 sometimes incorrectly 
interpreted published records. Moreover they used the original 
uncorrected record of HS98, while numerous corrections have 


been made to that during the last 17 years (e. 

g. Vaquero et al., 

2011; Vaquero & Trigo, 2014; Carrasco et al. 

, 2015b. Here we 


discuss some of the errors in ZP15, as examples of erroneous 
interpretation of historical data, in detail. 


2.4.1. Sunspot drawings vs. textual notes 


According to ZP15 “sunspot drawings provide a significantly 
larger number of sunspots, compared to textual or tabular 
sources”. This is trivial considering that the tabular sources 
often are related to astrometric observations of the Sun such 
as solar meridian altitude or the apparent solar diameter 
(IVaquero & Gallegol 120141) . However, if one considers only 
those tabular sources that contain explicit information about the 
presence or absence of sunspots then drawing s ources appea r to 
be consistent with the reliab le tabular sources dKovaltsov et all 
2004ICarrasco et all 12015b . 

The main assumption in ZP15 is that sunspots were omit¬ 
ted, especially in verbal reports, if they were not round and did 
not resemble a plane t. The only direct example o f that is given, 
with a reference to IVaquero & Vazquez! ( 2009 ), that Harriot 
drew three sunspots on Dec 8 (= Dec 18 greg.) 1610 but wrote 
that the Sun was “clear”. However, this discussion was based 
on an incorrect interpretation of the original texts. The actual 
statement of lHarriotl (1613) is 


The altitude of the sonne being 7 or 8 degrees. It being 
a frost and a mist. I saw the sonne in this manner [draw¬ 
ing]. I saw it twise or thrise. Once with the right ey and 
other time with the left. In the space of a minute time, 
after the sonne was to cleare. 


As indicated by the observing times and numerous other state¬ 
ments about a “well tempered” Sun in the course of his ob¬ 
servations, he mostly observed near sunrise or sunset or with a 
certain cloud cover to be able to look through the telescope. The 
statement that “the sonne was to cleare” refers to the fact that 
the Sun became too bright after a few minutes of observing. In 
this context, “cleare” means “bright” and not clear or spotless. 
Therefore, this example was incorrectly taken by ZP15 as an il¬ 
lustrative case of the discrepancy between textual and drawing 
sources. 

As another example, we compared the textual records 
by ISmogulecz & Schonberger ( 162 6), who had conservative 
views on sunspots (see Sect. 12.2. H) , with the drawings made by 
Schemer in Rome for the same period of 1625. We found that 
that Smogulecz and Schonberger omitted a number of spots 
from the drawings, but mentioned all spots they saw in their text 
(calling them ‘stellae’), in accordance with Scheiner. This is in 
contradiction with the assumption of ZP15 that verbal reports 
are subject to withholding spots. Table Q] lists the numbers of 
spots mentioned in their text versus those drawn in the figures. 
(We note that the value^are also incorrectly used in HS98.) 
Smogulecz & Schonberger ( 162 6) selected certain spots which 
were visible long enough to measure the obliquity of the Sun’s 
axis with the ecliptic and plotted them schematically as circles 
as not particularly interested in their shape. 


2.4.2. Relation between maximum number of sunspot 
groups and sunspot number 

ZP15 proposed a new method to assess the amplitude of the so¬ 
lar cycle during the MM. As the amplitude of a sunspot cycle, 
A * g , they used the maximum daily number of sunspot groups 
G* during the cycle, so that A G = 12.08 x G*, where the coeffi¬ 
cient 12.08 is a scaling between the average number of_sunspot 
groups and the sunspot number dHovt & Schattenl 1 19981) . We 
note that using the maximum daily value of G* instead of the 
average value G leads to a large overestimate of the sunspot 
cycle amplitude, particularly during the MM. We analyze the 
HS98 database for the period 1886-1945, when sunspot cycles 
were not very high, to compare the annually averaged group 
sunspot numbers Rg and the annual values of A G obtained us¬ 
ing the annual maxima of the daily sunspot group numbers G*. 
Fig. [5^ shows a scatter plot of the annual values of Rg and 
G* (dots), while the dashed line gives an estimate of the A* g 
based on G*, following the recipe of ZP15. One can see that 
while there is a relation between annual Rg and G*, the pro¬ 
posed method heavily overestimates the annual sunspot activ¬ 
ity. Fig. [5j> shows the overestimate factor Y = A* g /Rg of the 
sunspot numbers as a function of G*. While the factor Y is 2-3 
for very active years with G* > 15, the overestimate can reach 
an order of magnitude for years with weaker activity, such as 
during the MM. When applying to the sunspot cycle amplitude, 
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Table 1. Comparison of the number of spots listed in the ver¬ 
bal reports versus^ the number of spots in the drawings by 
Smogulecz & SchonbergeJ ( 16261) in 1625. 


Date 

Text 

Drawn 

Date 

Text 

Drawn 

1625 Jan 14 

1 

1 

1625 Aug 22 

2 

1 

1625 Jan 15 

1 

1 

1625 Aug 23 

2 

1 

1625 Jan 16 

4 

1 

1625 Aug 27 

6 

1 

1625 Jan 17 

8 

1 

1625 Aug 28 

10 

1 

1625 Jan 18 

2 

1 

1625 Aug 31 

7 

1 

1625 Jan 19 

4 

1 

1625 Sep 01 

6 

1 

1625 Jan 20 

2 

1 

1625 Sep 05 

8 

4 

1625 Feb 12 

8 

1 

1625 Sep 07 

6 

3 

1625 Feb 16 

10 

1 

1625 Sep 08 

6 

3 

1625 Feb 17 

11 

1 

1625 Sep 11 

5 

3 

1625 Feb 18 

10 

1 

1625 Sep 12 

4 

3 

1625 Feb 21 

4 

1 

1625 Sep 13 

2 

2 

1625 Jun 01 

9 

1 

1625 Oct 05 

9 

8 

1625 Jun 04 

3 

1 

1625 Oct 06 

2 

1 

1625 Jun 05 

3 

1 

1625 Oct 09 

4 

4 

1625 Jun 06 

2 

1 

1625 Oct 10 

7 

8 

1625 Jun 07 

3 

1 

1625 Oct 11 

9 

9 

1625 Jun 09 

2 

1 

1625 Oct 13 

2 

1 

1625 Aug 08 

6 

1 

1625 Oct 14 

2 

1 

1625 Aug 09 

4 

1 

1625 Oct 15 

3 

1 

1625 Aug 10 

2 

1 

1625 Oct 25 

1 

1 

1625 Aug 12 

4 

2 

1625 Oct 26 

1 

1 

1625 Aug 13 

3 

2 

1625 Oct 27 

1 

1 

1625 Aug 14 

3 

2 

1625 Oct 28 

1 

1 

1626 Aug 15 

4 

2 

1625 Oct 29 

1 

1 

1625 Aug 17 

2 

2 

1625 Oct 31 

1 

1 

1625 Aug 18 

4 

2 

1625 Nov 01 

1 

1 

1625 Aug 19 

2 

1 





the error becomes even more severe. Thus, by taking the cycle 
maxima of the daily number of sunspot groups instead of their 
annual means, ZP15 systematically overestimated the sunspot 
numbers during the MM by a factor of 5-15. 

The number of sunspot groups in 1642. 

ZP15 proposed that the solar cycle just before the MM was 
high (sunspot number « 100) which is based on a report of 
8 sunspot groups observed by Antonius Maria Schyrleus of 
Rheita in February 164 2 as presented in the HS 98 database. 
However, as shown bv iGomez & Vaquerol ( 2015b . this record 
is erroneous in the HS98 database because it is based on an in¬ 
correct translation from the original Latin records, which says 
that one (or a few) group was observed for 8 days in June 1642 
instead of 8 groups in February 1642. Accordingly, the max¬ 
imum daily number of sunspot groups reported for that cycle 
G* was 5, not 8, reducing the cycle amplitude claimed by ZP15 
(see Sect. 12.4.21) by about 40%. 

The number of sunspot groups in 1652. 

The original HS98 record contains 5 sunspot groups for the 
day of Apr-01-1652, referring to observations by Johannes 
Hevelius. Accordingly, this value (the highest daily G* for the 
decade) was adopted by ZP15 leading to the high proposed 
sunspot cycle during t he 1650s. However, as discussed by 


Vaquero & Trigol d2014l) in great detail, this value of 5 sunspot 


groups is an erroneous interpretation, by HS98 with reference 



Maximum number of groups G* 

Fig. 5. Illustration of the incorrectness of the method used by 
ZP15 to assess the group sunspot number Rq during the MM. 
Panel (A): Annual values of Rg as a function of the maximum 
daily number of sunspot groups G* for the same year in the 
HS98 catalogue for the period of 1886-1945; the dashed line 
is the dependence of A* G = 12.08 x G* used by ZP15. Panel 
(B): The overestimate factor Y of the GSN by the ZP15 method 
Y = A* g /R g . 


to [Wolf! ( 1856 ), of the original Latin text by Hevelius say¬ 
ing that there were 5 spots in two distinct groups on the Sun. 
Accordingly, the correct value of G* for that day should be 2 
not 5. 

The number of sunspot groups in 1705. 

A high sunspot number of above 70 was proposed by ZP15 
for the year 1705 based on six sunspot groups reported by J. 
Plantade from Montpellier (the correction factor for this ob¬ 
server is 1.107 according to HS98) for the day Feb-13-1705. 
This observer was quite active with regular observations dur¬ 
ing that period, with 44 known daily observational reports for 
the year 1705. For example, J. Plantade reported 2, 3, 6, and 
1 groups, respectively, for the days of Feb-11 through Feb-14. 
His reports also mention the explicit absence of spots from the 
Sun after the group he had followed passed beyond the limb. 
However, he did not make any reports during long spotless peri¬ 
ods, and wrote notes again when a new sunspot group appeared. 
The average number of sunspot groups per day reported by J. 
Plantade for 1705 was 1.22, which is a factor of 5 lower than 
that adopted by ZP15 who only took the largest daily value (see 
Sect. 12.4.21) . If one calculates the group sunspot number from 
the dataset of J. Plantade records for 1705 in the classical way, 
one obtains a value of Rg = 16.3 (1.22 x 12.08 x 1.107). 


2.5. Butterfly diagram 

According to sunspot drawings during some periods of 
the MM, a hint of the butterfly diagram has been iden¬ 
tified, particularly towards the end of the MM after 
1670 dRihes & Nesme-Rihesl 1 19931: ISoon & Yaskdll 12003 


Casas et al., 2006b. However, the latitudinal extent of the but 


terfly wings was quite narrow, being within 15° for the core 
MM (1645-1700) and 20° for the period around 1705, while 
cycles before and after the MM had a latitudinal extent of 30° 
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or greater. This suggests that the sunspot occurrence during the 
MM was limited to a more narrow band than outside the MM. 
Here we compare the statistics of the latitudinal extend of the 
butterfly diagram wings for solar cycles # 0 through 22 (cy¬ 
cles 5 and 6 are missing). The cycles #0 through 4 were cov¬ 
ered by digitized drawings made by Staudacher for the period 
1749-1792 (fXTitl 200 8 ). cycles #7 through 10 (1825-1867) 
were covered by digitized drawings made by Samuel Heinrich 


Schwabe (Arh et aL, 2013), cycle # 8 by drawings of Gustav 
Sporer foiercke et all 12015 ). while the period after 1874 was 


studied using the Royal Greenwich Observatory (RGO) cat¬ 
alogue. Moreover, a machine-readable version of the sunspot 
catalogues of the 19th century complied by Carrington, Peters 
and d e la Rue has been released recently dCasas & Vaquerol 
120141) . For each solar cycle we defined the maximum latitude 
(in absolute values without differentiating north and south) 
of sunspot occurrence. Since the telescopic instruments were 
poorer during the MM than nowadays, for consistency we con¬ 
sidered only large spots with the projected spot area greater 
than 100 msd (millionths of the solar disc). The result is shown 
in Fig. [6] as a function of the cycle maximum (in Rq). One can 
see that there is a weak dependence for stron ger cycles gen- 
erally having a larg er latitudinal span (cf, e.g . Ivitinskv et al 


erallv having a larger latitudinal span (cf, e.g . | 
1986; ISolanki et all 120081:Lliang et all 12011 hh . 
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but the latitu¬ 
dinal extent of the butterfly wing was always greater than 28° 
for the last 250 years. A robust link between the mean/range 
latitude of sunspot occurrence and cycle strength is related to 
the dynamo wave in the s olar convection zone an d has been 
empiric ally studied, e.g., bv ISolanki et al.1 d2008l) or l Jiang et al 


0 50 100 150 200 

Cycle maximum 

Fig. 6. Maximum latitudinal span of the butterfly diagram as 
a function of the cycle amplitude in annual Rg for solar cy¬ 
cles # 0-4 and 7-22, accounting only for large spots (with area 
greater than 100 msd). The dashed and dotted lines depict the 
maximum latitudinal extent of sunspot occurrence during the 
core (1645-1700) and the entire MM (1645-1712). 


be interpreted as sunspots have sometimes been used as 
an argument su g gesting high solar a ctivity during t he MM 


( Schovel 1983 : Nagovitsyn! 2001 : Ogurtsov et al. . 2003 


IZolotova & Ponvavinl 2015 ). Such statements are based on 


d2011al) . Since the maximum latitudinal extent of sunspots dur¬ 
ing the MM was 15° (during the core MM) or 20° (around ca. 
1705), it suggests a weak toroidal field causing a narrower lat¬ 
itudinal range of sunspot formation during the MM. This con¬ 
clusion is in ag reement with the result s of a more sophisticated 
analysis by llvanov & Miletskvl ( 2011 ) who found that the lat¬ 
itudinal span of the butterfly diagram during the late part of 
the MM should be 15-20°, i.e., significantly lower than during 
the normal cycles. One may assume that all the higher latitude 
spots were deliberately omitted by all the observers during the 
MM but we are not aware of such a bias. 

We note that two data sets of sunspot latitudes during the 
MM have been recently recovered and translated into machine- 
readable format (IVaquero et all l2015bh . Using these data sets, 
a decadal hemispheric asymmetry index has been calculated 
confirming a very strong hemispherical asymmetry (sunspots 
appeared mostly in the southern hemisphere) in the MM, as re¬ 
ported in earlier works dSporer , 1889; Ribes & Nesme-Ribesl 


an assumption that sunspots must be large to be observed 
and that this is possible only at a high level of solar activity. 
However, as shown below, this is not correct. While such his¬ 
torical records can be useful in a long-term perspective showing 
qualitatively the pre sence of several Grand minima during the 
last two millennia (Clark & Stephensoni 1 19781: IVaquero et all 
20021: IVaquero & Va zquez. 2009) including also the MM, this 


Il993l: IlSokololf & Nesme-Ribesl. 19941) . Another moderately 
asymmetric pattern wa s observed only in th e beginning of the 
Dalton minimum dArlil2009l:IUsoskin et alll2009bl) . Thus, this 
indicates that the MM was also a special period with respect to 
the distribution of sunspot latitudes. 


dataset is not useful for establishing the quantitative level of so¬ 
lar activity over short timespans due to the small number of in¬ 
dividual observations and/or the specific meteorological, soci¬ 
ological and historical conditions require d for such record (see 
Chapter 2 in lVaquero & VazquezL l2QQ9h . Moreover, it is very 
important to indicate that the quality of the historical record 
of naked-eye sunspot (NES) observations was not uniform 
through the ages (i.e. during the approximately two milliennia 
covered by the record). In fact, the quality of such records for 
the last four centuries was much poorer than that for the 12th- 
15th centuries, due to a change in the type of historical sources. 
In particular, the da ta coverage was reduced gr eatly after 1600 
(see Figure 2.18 in lVaquero & VazquezL 12009b . There are very 
few NES records during the century between the MM and the 
Dalton Minimum, representing the social conditions support¬ 
ing such observations and the maintenance of such records 
rather than sunspot activity. Therefore, the historical record of 
NES observations is not useful to estimate the level of solar 
activity during rec ent centuries (Eddy! 1 19831: iMossmanl 1 19891: 
Willis et allIl996h . 


2.6. East-Asian naked-eye sunspot observations 

East-Asian chronicles reporting observations for about two 
millennia, by unaided naked eyes, of phenomena that may 


2.6.1. Do NES observations imply high activity? 

It is typical to believe that historical records of NES observa¬ 
tions necessarily imply very high levels of solar activity (e.g., 
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Fig. 7. Probability density function for occurrence of the annual 
group sunspot numbers for the years 1848-1918. Panel A: The 
red solid line represents the years (50 years) without naked- 
eye spot (no-NES) reports, while the blue dotted line represents 
only the years (21 years) with NES. Panel B: The difference 
between the no-NES and NES probability density functions. 
Error bars represent the lcr statistical uncertainties. 


naked-eye sunspot observ ations during the period 1874-1918 


Ogurtsov et all 12003 ). assuming that observable spots must 
have a large area exceeding 1900 msd (millionths of the solar 
disc) with a reference to lWittmannl |l978). However, the latter 
work does not provide any argumentation for such a value and, 
as shown below, this is not correct. 

Here we show that reports of NES observations do 
not necessarily correspond to high activity or even to big 
spots. We compared the East-Asian sunspot catalogue by 
Yau & Stephensonl (ll988l) for the period 1848-1918 (25 re¬ 
ported naked-eye observations during 21 years) with data from 
the HS98 catalogue. Figure |7 ] shows the probability density 
functions (pdf’s) of the sunspot numbers for the years with and 
without NES observations. One can see that the probability of 
NES reports to occur does not depend on the actual sunspot 
number as the blue dotted curve in panel A is almost flat, while 
intuitively it should be expected to yield larger probability for 
high sunspot numbers and to vanish for small sunspot numbers. 
Moreover, there is no statistically significant difference in the 
sunspot numbers between the two pdfs. Accordingly, the null 
hypothesis that the two pdfs belong to the same population can¬ 
not be statistically rejected. Obviously, there is no preference 
to NES observations during the years of high sunspot numbers. 
The naked-eye reports appear to be distributed randomly, with¬ 
out any relation to the actual sunspot activity. Accordingly, the 
years with unaided naked-eye sunspot reports provide no pref¬ 
erence for the higher sunspot number. 

Next we study the correspondence between the NES 
reports and actual sunspots during the exact dates of 
NES observations (allowing for 1 day dating mismatch 
because of the local time conversion). The data 


dYau & Stephensonl 1 19881) . The dashed and dotted lines depict 


the 100% (all spots above this line are visible) and 0% (de¬ 
tectability threshold) probability of observin g a s unspot of the 
given area by an unaided ey e, according to lSchaefe3 ( 19931) and 
IVaquero & Vazque3 (120091) . 


the sunspot area were taken from the Royal Greenwich 
Observatory (RGO) sunspot group photographic catalogue 
(http://solarscience.msfc.nasa.gov/greenwch.shtml). Figure [8] 
shows, as filled circles, the largest observed sunspot area 
for the days when East-Asian NES observation were re¬ 
ported during the years 1874-1918 dYau & Stephe nson. 


Detectability limits of the NES observations (Schaefer 


1993 


Vaquero & Vazque3. 120091) are shown as dotted (no spots 


on 


smaller than « 425 msd can be observed by the unaided eye) 
and dashed (all spots greater than « 1240 msd are observable) 
lines. One can see that half of the reported NESs lie below or at 
the lower detectability limit and are not visible by a normal un¬ 
aided human eye, likely being spurious or misidentified records 
(cf. lWillis et idll 19961) . 

As an example we consider two dates with NES records 
with the smallest sunspots. A sunspot was reported to be seen 
by naked eyes on Feb-15-1900, when there were no sunspots 
on the Sun according to RGO, while there was one very small 
group (11 msd area) on the pervious day of Feb-14-1900. 
Another example of a NES report is for the day of Jan-30-1911 
when there was a single small group (area 13 msd) on the Sun 
(see also Fig. 9 in i Yau & Stephensonl 119881) . Such small groups 
cannot be observed by an unaided eye. Moreover, in agreement 
with the above discussion, even for big spots above the 100% 
detectability level, the relation to solar activity is unclear. Open 
dots in Figure [8] denote the area of the largest spot observed 
each year vs. the mean annual sunspot number for years 1874- 
2013. One can see that the occurrence of a large sunspot de¬ 
tectable by naked eye does not necessarily correspond to a high 
annual sunspot number, as it can occur at any level of solar ac¬ 
tivity from Rg = 3 to 200. 
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Annual number of sunspot groups 

Fig. 9. Probability density function of the occurrence of NES 
records by S.H. Schwabe as a function of the annual groups 
sunspot number during the years of Schwabe’s observations. 


Another data set is provided by the naked-eye observa¬ 
tions by Samuel Heinrich Schwabe, who recorded telescopic 
sunspot data in 1825-1867, but also occasionally reported on 
naked-eye visibilities of sunspot groups (Pavai et al., 2015, in 
prep.). We analyzed the (annual) group sunspot numbers for 
each event when Schwabe reported a naked-eye visibility, as 
shown in Fig. |9] in the form of a probability density function 
versus the annual group sunspot number. The NES reports were 
quite frequent during years with low sunspot activity (~ 25% of 
naked eye spots were reported for the years with Rg below 20, 
some of them even below 10). It is interesting that about 20% 
of naked eye observations by Schwabe were reported on days 
with a single group on the solar disc. We note that Schwabe 
certainly looked without the telescope when he saw a big group 
with it, so the selection may be biased towards larger spots. On 
the other hand, it is unlikely that he would watch out for naked- 
eye spots only if there were just few group on the Sun, therefore 
we do not expect an observational bias towards lower activity 
periods. 

Thus, a significant part of the East-Asian NES observa¬ 
tional reports are unlikely to be real observations and, even 
if they were correct, they do not imply a high level of solar 
activity. This implies that the NES reports cannot be used as 
an index of sunspot activity in a simple way (cf.jEddvl Il983 ; 


UIl AAAVAW./V VI uwiliupv/t WVtl Y AAA U. JAAAA^AV VV U.J yVA . j J-/V 

Willis et all 1 1 9961: TMossmanlll 9891: lusoskinl 12013 ). 


There are several NES records also during the MM but they 
are more rare (8 years during 1645-1715), as discussed in detail 
here. 

Three NES observations are reported for the years 1647, 
1648 and 1650, respectively, which fall in a long gap (1646- 
1651) of telescopic observations where only a generic state¬ 
ment by Hevelius exists. The exact level of sunspot activity 
during these years is therefore unknown. 

A NES report dated Apr 30 (greg.) 1655 falls in a small gap 
in the HS98 database but there is some activity reported in the 
previous month of March. The mean annual GSN for the year 
1655, estimated in the ‘loose’ model of IVaquero et al.1 d2015al) . 
is Rq = 5.7. 

A NES was reported in the Spring of 1656 which overlaps 
with a sunspot group reported by Bose in February. R g is esti¬ 
mated in the ‘strict’ model as 12.7 dVaquero et all 2015 ah . 

There are four NES records for the year 1665, but three of 
them are likely to be related to the same event in late February, 
and one to Aug 27, thus yielding two different observations. 
These events again fall into gaps in direct telescopic observa¬ 
tions with only generic statements available. For this year, only 
nine daily direct telescopic records, evenly spread over the sec¬ 
ond half-year, exist. The observers Hevelius and Mezzavacca, 
both claimed the absence of spots. The exact level of activity 
for this year is therefore unknown. Probably, there was some 
activity in 1665 but not high, owing to the direct no-spots 
records (cf. Sect. 12.2.2b . 

Another NES was reported to be observed for three days in 
mid-March 1684, which falls on no-spot records by la Hire. We 
note that this year was well covered by telescopic observation, 
especially in the middle and late year, and it was relatively ac¬ 
tive Rq = 11.7 dVaquero et all 12015a. ‘strict’ model). Account 
for the probably missed spot in March would raise the annual 
GSN value of this otherwise well observed year by less than 2. 

One more NES record is for the year 1709 (no date or 
even season given). That year was well observed by differ¬ 
ent observers, with some weak activity reported intermittently 
throughout the year. The mean annual GSN in the ‘strict’ model 
is Rq = 5.3. 

Thus, except for the year 1684, there is no direct clash be¬ 
tween the East-Asian NES records and European telescopic ob¬ 
servations, and the former do not undermine the low level of 
solar activity suggested by the latter. 


2.6.2. NES observations around the MM 

According to the well- established catalogue of NESs 


( Yau & Stephenson! Il988l) from Oriental chronicles, NESs 
were observed relatively frequently before the MM - 16 years 
during the period 1611-1645 are marked with the NES records. 
A direct comparison between the NES catalogue and the HS98 
database (with the correction by IVaquero et all 1201 lb shows 
that the NES records either are confirmed by Eurpean tele¬ 
scopic observations (Malapert, Schenier, Mogling, Gassendi, 
Hevelius) or fall in data gaps (after removing generic state¬ 
ments from the HS98 database). There is no direct contradic¬ 
tions between the datasets for that period. 


3. Indirect proxy data 

3.1. Aurorae borealis 

3.1.1. Geomagnetic Observations 

In recent years we have learned a great deal from geomag¬ 
netic observations about centennial-scale solar variability and 
how it influences the inner heliosphere, a nd hence the Earth 
( Lockwood et all 1 19991: iLockwoodi 12013 ). Such studies can¬ 
not tell us directly about the MM because geomagnetic activity 
was first observed in 1722 by George Graham in London and 
the first properly-calibrated magnetometer was not introduced 
until 1832 (by Gauss in Gottingen). Graham noted both regu¬ 
lar diurnal variations and irregular changes during the peak of 
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solar cycle # -3 (ca. 1720), which was the first significant cycle 
after the MM. This raises an interesting question: were these 
observations made possible by Graham’s advances to the com¬ 
pass needle bearing and observation technique or had magnetic 
activity not been seen before because it had not been strong 
enough? However, despite coming too late to have direct bear¬ 
ing on understanding the MM, the historic geomagnetic data 
have been extremely important because they have allowed us to 
understand and confirm the link between sunspot numbers and 
cosmogenic isotope data. In particular, they have allowed mod¬ 
elling of the open solar flux which shows that the low sunspot 
numbers in the MM are quantitatively (and not just qualita¬ 
tively) consistent with the high cosmogenic isotope abundances 


( Solanki et all 12000:10 wens et all 120121: Lockwood & Owena. 

20141) . This understanding has allowed the analysis presented 
in section [3741 


3.1.2. Surveys of historic aurorae 


Earlier in the same solar cycle as Graham’s first geomagnetic 
activity observations, on the night of Tuesday 17th March 1716 
(Gregorian calendar: note the original paper gives the Julian 
date in use of the time which was 6th March), auroral displays 
were seen across much of northern Europe, famously reported 
by Edmund lHallevI ( 1716 ) in Great Britain. 


What is significant about this event is that v ery few people 
in the country had seen an aurora before (IFaraLll996f) . Indeed, 
Halley’s paper was commissioned by the Royal Society for this 
very reason. This event was so rare it provoked a similar re¬ 
view under the auspices of l’Academie des Sciences of Paris 
(by Giacomo Filippo Maraldi, also known as Jacques Philippe 
Maraldi) and generated interest at the Royal Prussian Academy 
of Sciences in Berlin (by Gottfried Wilhelm Leibnitz). All these 
reviews found evidence of prior aurorae, but none in the previ¬ 
ous half century. 


Halley himself had observed the 1716 event (and correctly 
noted that the auroral forms were aligned by the magnetic field) 
but had never before witnessed the phenomenon. It is worth 
examining his actual words: “...[of] all the several sorts of me¬ 
teors [atmospheric phenomena] I have hitherto heard or read 
of, this [aurora] was the only one I had not as yet seen, and of 
which I began to despair, since it is certain it hath not happen’d 
to any remarkable degree in this part of England since I was 
born [1656]; nor is the like recorded in the English Annals since 
the Year of our Lord 1574.” This is significant because Halley 
was an observer of astronomical and atmospheric phenomena 
who even had an observatory constructed in the roof of his 
house in New College Lane, Oxford where he lived from 1703 
onwards. In his paper to the Royal Society, Halley lists reports 
of the phenomenon, both from the UK and abroad, in the years 
1560, 1564, 1575, 1580, 1581 (many of which were reported 
by Brahe in Denmark), 1607 (reported in detail by Kepler in 
Prague) and 1621 (reported by Galileo in Venice and Gassendi 
in Aix, France). Strikingly, thereafter Halley found no cred¬ 
ible reports until 1707 (Rpmer in Copenhagen and Maria and 
Gottfried Kirch in Berlin) and 1708 (Neve in Ireland). He states 
“And since then [1621] for above 80 years, we have no account 



Fig. 10. (a) The grey histogram shows the number of auroral 



nomena (including auroras) was only recently discovered and 
was compiled completely independently of other catalogues. 
Yet it shows, like the others, the dearth of sightings during the 
Maunder minimum, some events in 1707 and 1708 and the re¬ 
turn of regular sightings in 1716. (b) Annual variation of Na in 
the same dataset and of Ro¬ 


of any such sight either from home or abroad”. This analysis 
did omit some isolated sightings in 1661 from London (re¬ 
ported in the Leipzig University theses by Starck and Friiauff). 
In addition to being the major finding of the reviews by Halley, 
Miraldi and others (in England, France and Germany), a simi¬ 
lar re-appearance of aurorae was reported in 1716-1720 in Italy 
and in New England dSiscoelll980l) . 


The absence of auroral sightings in Great Britain during the 
MM is even more extraordinary when one considers the effects 
of the secular change in the geomagnetic field. For example, us¬ 
ing a spline of the IGRF (International Geomagnetic Reference 
Field, http://www.ngdc.noaa.gov/IAGA/ vmod/igr f.html) 
model after 1900 with the gufml model djackson et all 2000 ) 
before 1900 we find the geomagnetic latitude of Halley’s 
observatory in Oxford was 60.7° in 1703 and Edinburgh was at 
63.4°. Auroral occurrence statistics were taken in Great Britain 
between 1952 and 1975, and of these years the lowest annual 
mean sunspot number was 4.4 in 1954. Even during this low 
solar activity year there were 169 auroral nights observed at 
the magnetic latitude that Edinburgh had during the MM and 
139 at the magnetic latitude that Oxford had during the MM 
(Paton, 1959 ). In other words, The British Isles were at the 
ideal latitudes for observing aurora during the MM and yet 
the number reported was zero. This is despite some careful 
and methodical observations revealed by the notebooks of 
several scientists: for example, Halley’s notebooks regularly 
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and repeatedly use the term “clear skies” which make it 
inconceivable that he would not have noted an aurora had it 
been present. 

Halley’s failure to find auroral sightings in the decades be¬ 
fore 1716 is far from unique. Figure ITOl is a plot of auroral oc¬ 
currence in Great Britain from a previously unknown source. It 
is shown with the group sunspot number Rg during the MM. 
This catalogue of auroral sightings in the UK was published 
in 1870 by an astronomer and a Fellow of the Royal Society, 
E. J. Lowe, who used parish records, newspaper reports and ob¬ 
servations by several regular observers. His personal copy of 
the book (with some valuable ’’marginalia” - additional notes 
written in the margin) was recently discovered in the archives 
of the Museum of English Rural Life at the University of 
Reading, UK dLowe! Il870h . Here we refer to this personally 
commented copy of the book. We here have added to Lowe’s 
catalogue of English recordings the observations listed in the 
diary of Thomas Hughes in Stroud, England (Harrison, 2005) 
and the observations made by John Dalton in Kendall and later 
Manchester dDaltonlfl834l) . Like so many other such records, 
the time series of the number of auroral nights during each year, 
shown by the grey histogram in Figure [TOL reveals a complete 
dearth of auroral sightings during the MM. As such, this record 
tells us little that is not known from other surveys; however, 
it is important to note that this compilation was made almost 
completely independently of, and using sources different from 
other catalogues such as those by de Mairan and Fritz (see be¬ 
low). 

Figure [TOb shows the annual variation of the number of 
auroral nights and reveals the semi-annual variation (ISabine! 
1852) (equinoctial peaks in auroral occurrence were noted by 
de Mairail 1 173 3b . A corresponding semi-annual v ariation in 
geomagnetic activity dSabinel 1 18521: ICortiel Il912 ) is mainly 
caused by the effect of solar illumination of the nightside auro¬ 
ral current el e ctro je ts deliver et all l2000l: iLvatskv et all l200ll: 
Newell et~all l2002h . leading to equinoctial maxima in geo¬ 
magnetic activity. On the other hand, lower-latitude aurorae 
are caused by the inner edge of the cross-tail current sheet 
being closer to the Earth, caused by larger open flux in the 
magnetosphere-ionosphere system (iLockwood! 12013b and so 
are more likely to be caused by the effect of Earth’s dipole 
tilt on solar wind-magnetosphere coupling and, in particular 
the magnetic reconnection in the magnetopause that generates 
the open flux IRussell & McPherronlll973u . This is convolved 
with a summer-winter asymmetry caused by the length of the 
annual variation in the dark interval in which sightings are pos¬ 
sible. Note that FigureflOb shows a complete lack of any annual 
variation in group sunspot number, as expected. This provides 
a good test of the objective nature of the combined dataset used 
in Figure [TO) Both parts (a) and (b) of Figure [T0| are very simi¬ 
lar in form to the corresponding plots made using all the other 
catalogues. 

Elsewhere, however, other obser vers in 1716 were famil - 
iar with the phenomenon of aurorae dBrekke & Egeland . 1994). 
For example Joachim Ramus, a Norwegian (born in Trondheim 
in 1685 but by then living in Copenhagen), also witnessed au¬ 
rora in March 1716, but unlike Halley was already familiar 
with the phenomenon. Suno Arnelius in Uppsala had written 


a scientific thesis on the phenomenon in 1708. Indeed after 
the 1707 event Rpmer had noted that, although very rare in his 
native Copenhagen, such events were usually seen every year 
in Iceland and northern N orway (although it is not known on 
what basis he stated this) dStauningl 1201 ll: iBrekke & Egeland . 


1994). But even at Nordic latitudes aurorae appear to have 


been relatively rare in th e second half of the 17th century 
dBrekke & Egeland! 1 19941) . Petter Dass, a cleric in Alstahaug, 
in middle Norway, who accurately and diligently reported all 
that he observed in the night sky between 1645 and his death 
in 1707, and had read many historic reports of aurorae, never 
once records seeing it himself. In his thesis on aurorae (com¬ 
pleted in 1738), Peter Mpller of Trondheim argues that the au¬ 
rora reported over Bergen on New Year’s Eve 1702 was the 
first that was ever recorded in the city. Celsius in Uppsala was 
15 years old at the time of the March 1716 event but subse¬ 
quently interviewed many older residents of central Sweden 
and none had ever seen an aurora before. Johann Anderson was 
the mayor of Hamburg and discussed aurorae with Icelandic 
sea captains. They told him that aurorae were seen before 1716, 
but much less frequently (reported in iHorrebow! 1 1752 ). An im¬ 
portant contribution to the collation of reliable auroral observa¬ 
tions was written in 1731 by Jean-Jacques d’Ortous de Mairan 
dde Mairan! 1733 ), with a second edition published in 1754. 
Both editions are very clear in that aurorae were rare for at 
least 70 y ears before their r eturn in 1716. Th e m ore thorough 


surveys bv lLoveringl d 1 860h . lFritzl (fl873ll 188 1|) andl 


1978) have all confirmed this conclusion (see Eddy 


ink(1964 

19761) 


3.1.3. Reports of aurorae during the Maunder 
minimum 

The above does not mean that auroral sighting completely 
ceased during the Maunder minimum, de Mairan’s original 
survey reported 60 occurrences of aurorae in the interval 
1645-1698. Many authors have noted that the solar cycle in 
auroral occurrence continued during the Maunder minimum 


Unkt 1978; Vit.inskilL M97alGleissherglM9T7HSchrode^h992l: 


Legra nd et al. . 1 1 992l) . One important factor that must be con- 
sidered in this context is the magnetic latitude of the observa¬ 
tions. It is entirely possible that aurora were always present, at 
some latitude and brightness, and that the main variable with 
increasing solar activity is the frequency of the equatorward 
excursions of brighter forms of aurorae. In very quiet times, 
the aurora would then form a thin, possibly fainter, band at 
very high latitudes, with greatly reduced chance of observa¬ 
tion. An important indication that this was indeed the case 
comes from a rare voyage into the Arctic during the MM by 
the ships Speedwell and Prosperous in the summer of 1676. 
This was an expedition approved by the then secretary to the 
British Admiralty, Samuel Pepys, to explore the north east pas¬ 
sage to Japan. Captained by John Wood, the ships visited north¬ 
ern Norway and Novaya Zemlya (an Arctic archipelago north 
of Russia), reaching a latitude of 75° 59’ (geographic) before 
the Speedwell ran aground. Captain Wood reported that auro¬ 
rae were only seen at the highest latitudes by the local seaman 
that he met. Furthermore, anecdotal evidence was supplied by 
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Fritz who quoted a book on Greenland fisheries that auro 
were sometimes seen in the high Arctic at this time. The pi 
sibility of aurora watching at very high latitudes was the m 
criticism of de Mairan’s catalogue by Ramus, claiming tha 
relied on negative results from expe ditions that were outs 
the observing season set by sunlight (Brekke & Egela nd[[l9^~ 
Stauningil201lh . 


3.1.4. Comparison of Aurorae during the Maunder a 
Dalton minima 

The debate about the reality of the drop in auroral occurrer 
during the Maunder minimum was ended when a decline v 
seen during the Dalton minimum (c. 1790-1830). This mi 
mum is seen in all the modern catalogues mentioned above and 
in others, such as that bv iNevanlinnal ( 1995 ) from Finnish ob¬ 
servatories, which can be calibrated against modern-day obser¬ 
vations (Nevanlinna_& Pulkkinen, 200lb . Many surveys show 
the MM to be deeper than the Dalton minimum in auroral oc¬ 
currence but not by a very great factor (e.g.. lSilvermanLfl992l) . 
However, given the likelihood that aurorae were largely re¬ 
stricted to a narrow band at very high latitudes during both min¬ 
ima, observations at such high latitudes become vital in estab¬ 
lishing the relative depths of these two minima. In this respect 
the survey by Vazquez et al. (2014) is particularly valuable as, 
in addition to assigning locati ons to every sighti ng, i t includes 
the hig h latitude catalogue s by Rubenson ( 18821) andlTromholt 
d 18981) as well as those of lSilvermanl (1992) andlFritz ( 1873 ). 


The quality control employed by Vazquez et al. (120141) means 
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that their survey extends back to only 1700 which implies that 
it covers 15 years before the events of 1716 and hence only the 
last solar cycle of the MM. 

Figure m is an analysis of the occurrence of aurorae be¬ 
tween 1700 and 1900. The green line is the number of auro¬ 
ral nights per year at geomagnetic latitudes_below 56° from 
a com b inatio n of th e cat alogues of Nevanlinna J 1995 ). iFritzl 
dl873l) . IFritzl dl88ll) and iLegrand & Simonl dl987 ). This se¬ 
quence clearly shows that aurorae at geomagnetic latitudes be¬ 
low 56° were indeed rarer in both the last cycle of the MM 
and the two cycles of the Dalton minimum (DM). However, the 
number of recorded auroral sightings was significantly greater 
during DM than that in the MM. The points in Figure [HI show 
the geomagnetic latitude and time of auroral sightings from the 
catalogue of I Vazquez et al.1 d2014l) (their Figure 9). Black dia¬ 
monds, red squares and red triangles are, respectively, for ob¬ 
serving sites in Europe and North Africa, North America, and 
Asia. Blue dashed lines mark the minimum latitude of auroral 
reports in the last solar cycle of the MM and in the two cycles of 
the DM. During the Dalton minimum many more aurorae were 
reported (symbols in the Figure) poleward of the 56° latitude. 

Considerably fewer arcs were reported at the end of the 
MM at these latitudes, despite the inclusion by 1 Vazquez et al.1 
( 2014 ) of two extra catalogues of such events for this period 
at auroral oval latitudes. Furthermore, the two dashed lines 
show the minimum latitude of events seen in these two minima: 
whereas events were recorded down to a magnetic latitude of 
45° in the DM, none were seen at the end of the MM below 55°, 


Fig. 11. Occurrence of auroral reports, 1700-1900. The green 
line is the number of auroral nights at geomagnetic latitudes be¬ 
low 5 6° from a combina t ion of several catalogues dNevanlinnai 
Il99.it iFritTt 1 18731 Il88lt ILegrand & Si J Il987h . The points 
show the geomagnetic latitude and time of auroral sightings 
from the catalogue of 1 Vazquez et al. ( 12014 ) (their Figure 9). 
Black diamonds, red squares and red triangles are, respectively, 
for observing sites in Europe and North Africa, North America, 
and Asia. Blue dashed lines mark the minimum latitude of au¬ 
roral reports in the last solar cycle of the Maunder minimum 
(MM) and in the two cycles of the Dalton Minimum (DM). 


consistent with the dearth of observations in central Europe at 
the time. Note that during the MM there are some observa¬ 
tions at magnetic latitudes near 27°, all from Korea (with one 
exception which is from America). They were reported to be 
observed in all directions, including the South, and to be red 


(Yauetal., 1995) which makes them unfavorable candidates 


for classical aurorae. By their features these could have been 
stable auroral red (SAR) arcs (Zhang, 1985) which in mod¬ 
ern times are seen at mid-latitudes mainly during the recovery 


phase of geomagnetic storms (Kozyra et al., 1997). These arcs 


are mainly driven by the ring current and differ from thejior- 
mal auroral phenomena. Moreover, as stated bv lzhangl ( 1985 ) 
“We cannot rule out the possibility that some of these Korean 
sightings were airglows and the zodiacal light”. We here con¬ 
centrate on the higher latitude auroral oval arcs. Note that the 
plot also shows the return of reliable lower latitude sightings in 
Europe in 1716 and in America in 1718. 

Figure [12] corresponds to Figure \T0\ but now based 
on a compilation of all major historical auroral catalogues. 
Figure fl2l employs the list of aurora days bv lKrivskv & Peiml 
( 1988 ) which is based on 39 different catalogues of obser¬ 
vations at geomagnetic latitudes below 55° in Europe, Asia, 
North Africa. To this has been added the catalogue of lLoveringi 
d 18671^1 of observations made in and around Cambridge, Mass., 
USA which was at a magnetic latitude close to 55° in 1900, and 
the recently-discovered catalogue of observations from Great 
Britain by Lowe ( 1870 ). Figure fl2h shows the low level and 
gradual decline in the occurrence of low- and mid-latitude au¬ 
roral observations during the MM and in the decades leading 


2 Paper and catalogue are available from 
http: //w w w.j stor. org/stable/25057995 
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Fig. 12. Same as Figure[lOl but compiled from 41 different cat¬ 
alogues of auroral observations at magnetic latitudes below 55° 
in Europe, Asia, North Africa, New England and Great Britain. 
The time series covers both the Maunder and the Dalton min- 


Rilev et all l2Q15h suggest that the sol ar corona w as reddish 


in 


Riley et al. (120151) 


3.3. Heliospheric conditions 

The Sun was not completely quiet during the MM, and a 
certain level of heliospheric activity was still present - the 
heliosphere existed, the solar wind was blowing, the helio¬ 
spheric magnetic field was there, although at a stro ngly re 
duced level (e.g. Cliveretail 19981: ICaballero-Lopez etal 


l2004 iMcCracken & Beeil 2014 ). Since heliospheric distur¬ 


bances, particularly those leading to cosmic ray modulation, 
are ultimately driven by solar surface magnetism (iPotgieter . 
2013) , which is also manifested through sunspot activity, cos¬ 


tma. 


up to it. This is in contrast to the general rise in observation 
reports that exists on longer timescales as scientific recording 
of natural phenomena became more common. After the MM, 
solar cycles in the auroral occurrence can clearly be seen and 
the correlation with the annual mean group sunspot numbers 
Rg is clear. Even for these lower latitude auroral observations 
it is unquestionable that the MM is considerably deeper than 
the later Dalton minimum. The annual variability (Figure fl2b) 
is obvious also for this data set. 

From all of the above it is clear that the Maunder minimum 
in auroral, and hence solar, activity was a considerably deeper 
minimum than the later Dalton minimum. 


3.2. Solar corona during the Maunder minimum 

As shown already bv lEddv ( 1976 ) and re-analyzed recently by 


mic ray variability is a good indicator of solar activity, espe¬ 
cially on time scales longer tha n a solar cycle (Beer, l200d: 
iBeer et all 1201 2l: lUsoskinl 120131) . Here we estimate the helio¬ 
spheric conditions evaluated for the period around the MM, 
using different scenarios of solar activity, and compare those 
with directly measured data on cosmogenic isotopes in terres¬ 
trial and extra-terrestrial archives. 

The open solar magnetic flux (OSF) is one of the main he¬ 
liospheric parameters defining the heliospheric modulation of 
cosmic rays. It is produced from surface magnetic fields ex¬ 
panding into the corona from where they are dragged out into 
the heliosphere by the solar wind. Consequently, it can be mod¬ 
elled using the surface distribution of sunspots and a model of 
the surface magnetic flux transport dWang & Sheelevl [2002). 
If their exact surface distribution is not known, the number 
of sunspots can also serve as a good input to OSF computa¬ 


tions, using models_ofmagnetic fluxevolution (ISolanki et al 


2000i 12002c iLockwood & Qwensll2014l) o r with more comple x 


surface flux transport simulations (e.g., Iliang et all 1201 Lb). 


Here we use the simpler, but nonetheless very successful 
model to calculate th e OSF from the sunspot number series 
( Lockwood & Owensl. 2014 : Lockwood et al. . 2014al) . This 
model quantifies the emergence of open flux from sunspot 
number using an analysis of the occurrence rate and mag¬ 
netic flux content of coronal mass ejections as a function of 
sunspot number over recent solar cycles dOwens & Lockwood . 


20121) . The open flux fractional loss rate is varied over the so- 


Riley et al. (12015 ). recorded observations of solar eclipses sug¬ 
gest the virtual absence of the bright structured solar corona 
during the MM. While 63 total solar eclipses should had taken 
place on Earth between 1645 and 1715, only four (in years 
1652, 1698, 1706 and 1708) were properly recorded in a sci¬ 
entific manner, others were either not visible in Europe or 
not described in sufficient detail. These reports (see details in 


and unstructured which was interpreted dEddv . 1976 ) as the F- 
corona (or zodiacal light) in the absence of the K-corona. The 
normally structured corona reappeared again between 1708 and 
1716, according to later solar eclipse observations, as discussed 


Observations of the solar corona during total eclipses, al¬ 
though rare and not easy to interpret, suggest that the corona 
was very quiet and had shrunk during the MM, with no large 
scale structures such as streamers. This also implies the decline 
of surface activity during the MM. 


lar cycle with the current sheet tilt, as predicted theoretically 
bv [Owens et al.l (201 If) and the start times of each solar cycle 
taken from sunspot numbers (during the MM the 10 Be cycles 
are used). The one free parameter needed to solve the continu¬ 
ity equation, and so model the OSF, is then obtained by fitting 
to the open flux r econstruction derived from geomagnetic data 
for 1845-2012 bv Lockwood et aD d2014al) . 

We computed OSF series (Figure [T3]) corresponding to the 
two scenarios for the number of sunspots during the MM (see 
Section O, viz. L- and H-scenarios. 

The OSF is the main driver of the heliospheric modula¬ 
tion of cosmic rays on time-scales of decades to centuries (e.g., 
lUsoskin et all 120021) . with additional variability defined by the 
heliospheric current sheet (HCS) t ilt and the large scale polarity 
of the heliospheric m agnetic field (lAlanko-Huotari et all 12006c 
homas et al. L 20141). Using a n updated semi-empirical m odel 


dAlanko-Huotari et all 120061: lAsvestari & Usoskim 12015k we 


have computed the modulation potential 0 (see definition and 
formalism in lUsoskin et al.l 12005 ) for the period since 1610 
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Fig. 13. Evolution with time of the open solar magnetic flux, 
OSF, (panel A) and of the modulation potential 0 (panel B) for 
the two scenarios of solar activity during the MM (see text). 


for the two scenarios described above, as shown in Figure [l3b. 
These series will be used for a subsequent analysis. 

3.4. Cosmogenic radionuclides 

The cosmogenic radionuclides are produced by cosmic rays in 
the atmosphere, and this forms the dominant source of these 
isotopes in the terrestrial system feeer et al 112012k . Production 
of the radionuclides is controlled by solar magnetic activity 
quantified in the heliospheric modulation potential (see above) 
and by the geomagnetic field, both affecting the flux of galactic 
cosmic rays impinging on Earth. For independently known pa¬ 
rameters of the geomagnetic field, one can use the measured 
abundance of cosmogenic radioisotopes in a datable archive 
to reconstruct, via proper modelling including production and 
transport of the isotopes in the Earth’s atmosphere, the so¬ 
lar/helio spheric magnetic activity in the past ( Beer et all 12012 : 
Usoskinl[2013 ). Here we used a rece nt archeomagneti c recon¬ 
struction of the geomagnetic field dEicht et all 12013 ) before 
1900. In the subsequent subsections we apply the solar modula¬ 
tion potential series obtained for the two scenarios (Figurefl3b) 
to cosmogenic radionuclides. 


Years 

Fig. 14. Time profile of decadally smoothed radiocarbon 14 C 
production rate. The black curve with grey error bars represents 
the reconstruc t ion by iRoth & Joosl d20 13 ) based on the Intcall 3 
dReimer et all 12013 ) radiocarbon data. Colored curves depict 
the computed production for the two scenarios as labeled in the 
legend (see text for details). 


In the same plot we show also the 14 C prod uction rate ob 
tained bv [Roth & Joosl ( 2013 ) from the Intcall3 dReimer et al 


3.4.1. 14 C in tree trunks 

Using the recent model of radiocarbon 14 C production 
( Kovaltsov et all 12012 ). we have computed the expected global 
mean radiocarbon production rate for the two scenarios ana¬ 
lyzed here, as shown in Fig.[l4l One can see that the variability 
of 14 C production is quite different in the H-(red dotted) and 
L-(blue solid curve) scenarios. While the former is rather con¬ 
stant, with only a weak maximum during the MM, even smaller 
than that for the Dalton minimum in the early 1800s, the latter 
exhibits a high and long increase during the MM, which is sig¬ 
nificantly greater than that for the DM in both amplitude and 
duration. 


2013k global radiocarbon data, using a new generation state- 
of-the-art carbon cycle model. The 14 C global production ex¬ 
pected for the L-scenario agrees very well with the data, within 
the uncertainties, over the entire period of 1610-1880 con¬ 
firming the validity of this scenario. On the contrary, the H- 
scenario both quantitatively and qualitatively disagrees with 
the observed production during the MM, implying that the so¬ 
lar modulation of cosmic rays is grossly overestimated during 
the MM by this scenario. Thus, the 14 C data support a very 
low level of heliospheric (and hence solar surface magnetic) 
activity during the MM, a level that is considerably lower than 
during the Dalton minimum. 


3.4.2. lu Be in polar ice cores 

With a similar approach to that taken for the analysis of 
14 C (section 13.4.11) we have computed the depositional flux 
of 10 Be in polar_ regions. We used the same archeomag- 
netic mod el dLicht et al. . 2013k , the r ecent 10 Be production 
model by [Kovaltsov & Usoskinl ( 2010 ). and the atmospheric 
transp ort/deposition model as parameterized bv lHeikkila et al.1 

(B . 

The results are shown in Figure |T5j As discussed in the 
previous section, the expected curve for the H-scenario (red 
dashed curve) shows little variability, being lower (implying 
higher solar modulation) during the MM than during the DM, 
while the L-scenario yields a higher flux (lower modulation) 
during the MM. The two grey curves depict 10 Be fluxes mea¬ 
sured in two opposite polar regions. One is the data series of 
10 Be depositional flux measured in the Antarctic Dome Fuji 
(DF) ice core dHoriuchi et all 112007 ). The other is the 10 Be flux 
series measured in the Greenland NGRIP (Nor th Greenland 


Ice-core Project) ice core dBerggren e t al., 200 9). Because of 


the different local climate conditions (Heik kila et al.ll200~9k . the 
latter was scaled by a factor of 0.83 to match the same level. 




















































































































Usoskin et al.: Maunder minimum: A reassessment 


17 



Fig. 15. Depositional flux of 10 Be in polar ice. All data are 
pseudo-decadal and ~ 25-yr smoothed because of strong high- 
frequency noise. Grey curves depict the fluxes measured at 
Dome Fuji (DF) in Antarctica ( Horiuc hi et al., 2007), and at 


NGRIP (scaled as 0.83), Greenland feerggren et all 2009 ). 


Error bars are estimates of both statistical and systematic er¬ 
rors. The blue and red curves depict the modelled 10 Be depo¬ 
sitional fluxes for the L- and H-scenarios (see text for details), 
respectively. 


This scaling does not affect the shape of the curve and in par¬ 
ticular not the ratio of the 10 Be flux in the MM and the DM. 
One can see that while the time profiles of the two datasets 
differ in detail, probably because of the different climate pat¬ 
terns iUsoskin et all 12009a ) and/or timing uncertainties, both 
yield high 10 Be production during the MM. This corresponds 
to extremely low solar activity dMcCracken et all 120041) . The 
L-scenario agrees with the data reasonably well (the data dis¬ 
play even higher maxima than the model), while the H-scenario 
clearly fails to reproduce the variability of 10 Be measured in 
polar ice. 

Thus, the 10 Be data from both Antarctic and Greenland ice 
cores support a very low level of heliospheric (and hence solar 
surface magnetic) activity during the MM, significantly lower 
than during the Dalton minimum. 


3.4.3. 44 Ti in meteorites 

While records of terrestrial cosmogenic radionuclides may 
be affected by transport an d deposition p r ocesses, which are 
not a lways exactly known ( Usoskin et all l2009al: iBeer et af 


20121) . cosmogenic nuclides measured in fallen meteorites are 
free of this uncertainty, since the nuclides are produced di¬ 
rectly in the meteorite’s body while in space, and measured 
after their fall on the Earth. However, time resolution is lost or 
at least greatly reduced in this case, and the measured activ¬ 
ity represents a balance between production and decay over the 
time before the fall of the meteorite. An ideal cosmogenic nu¬ 
clide for our purpose is 44 Ti with a half- life of about 60 years 


(Ahmad et al 


purpos( 

.,1998; 


Bon ino et all 1995). Here we test the cos¬ 


mic ray variability as inferred from different scenarios of solar 
activity since 1600, following exactly the method described in 
detail in IUsoskin et al . (200 6) and the dataset of 44 Ti activity 


measured in 19 stony meteorites fallen between 1776 and 2001 
(Taricco et al., 2006). Applying the modulation potential series 



Years 

Fig. 16. Time profile of the 44 Ti activity in units of disintegra¬ 
tions per minute per kg of iron and nickel in the meteorite. The 
grey dots with error bars are the measurements dTaricco et all 
2006b . The colored curves depict the computed activity for the 
two scenarios (see text). 


as described in Sect. 13.31 to the cosmic ray flux and calculat¬ 
ing the expected 44 Ti activity as a function of the time of the 
meteorite’s fall, we compare the model computations for the 
two scenarios with the measured values in Fig. Q2] One can 
see that the series for the L-scenario fits the data rather well, 
whereas the curve resulting from the H-scenario lies consider¬ 
ably too low. As a merit parameter of the agreement we use the 
X 1 value for the period of 3.3 half-lives (attenuation factor 10) 
since the middle of the MM, viz. until 1880, which includes 
7 meteorites fallen between 1776 and 1869 (7 degrees of free¬ 
dom). The^ 2 (7) for the L-scenario is 2 (0.33 per degree of free¬ 
dom), which perfectly fits the hypothesis. For the H-scenario, 
however, x 2 (7) = 17.6 (2.93 per degree of freedom) indicat¬ 
ing that this hypothesis should be rejected with a high signifi¬ 
cance of 0.014. Note that the long residence times of the me¬ 
teorites within the helisophere means that the difference during 
the MM between the two scenarios has an influence on the pre¬ 
dicted 44 Ti abundances even for meteorites that fell to Earth rel¬ 
atively recently. This results in the observed abundances being 
inconsistent with the H-scenario for a large number of the me¬ 
teorites whereas they are consistent, within the observational 
uncertainty, with the L-scenario. 

Accordingly, the hypothesis of a high level of solar activ¬ 
ity during the MM is rejected at a high significance level using 
indirect data of 44 Ti in meteorites, while the conventional sce¬ 
nario of very low activity during the MM is in full agreement 
with the data. 


4. Consequences of the Maunder minimum 

4.1. Solar/stellar dynamo 


Major changes in the secular level of solar activity, such as 
Grand minima/maxima form a challenge for our understand¬ 
ing of the origin and evolution of the solar magnetic field. 
It has been recently shown, by analyzing the sunspot num¬ 
bers reconstructed from 14 C for the last 3000 years, that Grand 
minima form a separate mode of solar activity ( Usoskin et all 
2014), likely corresponding to a special regime of the solar dy¬ 


namo. In addition to the traditional concept of cyclic solar ac- 
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tivity associated with periodic nonlinear oscillations of large- 
scale magnetic field, solar dynamo models include now, as 
their natural element, various deviations from pure periodicity. 
Accordingly, it is crucially important for our understanding of 
solar/stellar dynamos to know whether Grand minima indeed 
exist and what their parameters are. Although a possibility of 
direct modelling of the MM by dynamo theory is limited by the 
lack of information concerning flows in the solar interior, some 
important observational results pointing to quite peculiar fea¬ 
tures of the solar surface magnetic field configuration (slower, 
but more differential, rotation, strong hemispheric asymmetry 
of sunspot formation, and a possibly variable solar diameter) 


during the MM have been found (Ribe s & N esme-Ribes. 1993 
Sokoloff & Nesme-Ribes, 1994 ). 


Most large-scale solar dynamo models operate with aver¬ 
aged quantities taken as statistical ensembles of a moderate 
number of convective cells (see, e.g., a review by Charbonneau, 
l2010l) and include, in addition to the solar differential rota¬ 
tion, the collective inductive effect of mirror asymmetric con¬ 
vective turbulence and/or meridional circulation. The observa¬ 
tional knowledge of the turbulent quantities is especially lim¬ 
ited, they have t o be estimated using local direct numerica l 
simulations (e.g. ISchrinner et all 120051: iKapvla et~all 120091) . 
Generally, fluctuations of the intensity of the main drivers 
of the dynamo up to 10%-20% are expected. By includ¬ 
ing such fluct uations into a mea n -field dynamo mode l s it is 
possible (e.g. iMoss et all l2008bt IChoudhuri & Karaki 12012 : 
Passos et al. . 120141) to numerically reproduce a dynamo behav¬ 


forces and the momentum equation. In those cases, grand 
minima are produced w i thout a stochastic d ynamo effect 
dKiiker et all fl999l: IPipinl. 1 19991: iBushbvl 120061 for spherical 


shells). The presence of grand minima is therefore a natural fea¬ 
ture in mean-field dynamo modelling, including ‘side-effects’ 
like differential rotation variation and mixed parity, and is a re¬ 
sult of stochasticity and the nonlinearity of the full MHD equa¬ 
tions. 

Thus, the known phenomenology of the MM does not favor 
its interpretation as just a modulation of the normal 22-year 
Hale cycle by an additional longer cycle (Gleissberg cycle) as 
proposed, e.g., by ZP15. The very asymmetric and suppressed 
sunspot activity, accompanied by slower differential rotation, 
during the late phase of the MM implies, in the light of dynamo 
theory, a special mode of dynamo operation, leading to Grand 
minima. 

Meanwhile, it is still difficult to perform a direct nu¬ 
merical modelling of the MM, including sufficiently small 
scales to adequately model the convective turbulence (e.g. 
Charbonneaul 120101) . Although modern solar-type dynamo 
models reproduce the “regular” part of the behavior of 


the large-scale solar magnetic field dGhizaru et al 


2010 ; 


Brown et al. . l201 lHSchrinner et al.l.l2012l:lGastine et all 2012 : 


ior which deviates from the stable cyclic evolution and de¬ 
picts variability like the Grand minima. It is important that 
asymmetric magnetic configurations of mixed parity with re¬ 
spect to the solar equator can be e xcited even in the frame¬ 
work of conventional ctQ dynamo dBrandenburg et all 1 19891: 
Jennings & Weissl Il991|). similar to that existing on th e Sun 
at the end of the MM dSokoloff & Nesme-Ribes . 1994 ). The 


asymmetric sunspot occurrence during the late MM may be 
a signature of an unusual mode of the dynamo. It is known 
for spherical dynamos (Moss et all l2008aj) that even relatively 
moderate deviations from the nominal parameters associated 
with normal cycles can lead to the excitation of specific mag¬ 
netic configurations, for example, with a quadrupolar symme¬ 
try with respect to the solar equator. On the other hand, such an 
asymmetry is not expected for regular cycles with normal val¬ 
ues of the driving parameters. This, along with the suppression 
of the cycle amplitude, may be a specific feature of a MM-type 
event. 

An interesting fact is that the solar surface rotation was re¬ 
ported to be slower and more differential (changing faster with 
latitude) during the second half of the MM than during mod¬ 
ern times ( Ribes & Nesme-Ribesi, 1993 ). The enhanced dif¬ 
ferential rotation may have lasted until the mid-17th century 
dArlt & Frohlichll2012l) . These facts are related to the operation 
of the solar dynamo during the MM. Since the solar differen¬ 
tial rotation is a main driver of the dynamo, and the asymmetry 
implies a specific configuration, the dynamo could had being 
operating in a special state during this period. 

Differential rotation modulation and mixed parity are also 
an outcome of nonlinear dynamo models including Lorentz 


Kapvl a et al. . l2012h . it is still challenging to extend the integra¬ 
tion time over several magnetic cycles as needed to reproduce a 
Maunder-like minimum (see, however. [Passos & Charbonneau . 
2014 : lAugustson et al.[l2015 ). Thus, while the MM is identified 
as a special mode of the solar dynamo, we are not yet able to 
precisely model it. 

4.2. Solar irradiance 

It is now widely accepted that variations in solar irradiance in 
different wavelengths on time scales longer than about a day 
are driven by changes in the solar surface coverage by magnetic 
features, such as sunspots, that lead to a darkening of the solar 
disc, and facul ae or network e l ements, that lead to a brig hten- 
ing, (see, e.g.. iDomingo et all 120091: ISolanki et all l2013l and 
references therein). Consequently, the number of sunspots and 
faculae present during the MM would affect both the total (TSI) 
and the spectral (SSI) irradiance of the Sun at that special time 
in the history of solar activity. 

Direct observations of TSI/SSI are availabl e onl y 


from 1978 onwards (e.g., iFrohlichl 120131: iKoppl 120141) . 


Consequently, a number of models have been developed that 
reconstruct solar irradiance back to the MM. The TSI produced 
by a selection of such models is plotted in Fig.[T71 The selected 
models are based on very different data and techniques. 

Steinhilber et al.1 d2009h and iDelavgue & Bardl d201ll) ob¬ 
tain their reconstructions from timeseries of 10 Be concentra¬ 
tions in ice (green colour tones in the figure). Both studies used 
a simp le linear regression to calculate TSI. ISteinhilber et al. 
( 20098) first estimated the interplanetary magnetic field from the 


10 Be data and then used a lin e ar regression with the me asured 
TSI following |Frohlichld2009h . lDelavgue & Bardl d201ll) scaled 
the 10 Be record by assuming a fixed change in TSI between the 
MM and the last decades. Not shown in Fig. [IT] and also rely- 
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Fig. 17. Selected TSI re constructions sin c e 160 0, labeled 
in the plot are: Sea09 - ISteinhi lber et al. d2009b; DB11- 


Delaygue & Bard (2011); Wea05 


Krivovaetal. (20 


Veall - Vieira et al 


0); Deal4 


Wang et al. 


2005); KealO 
d2014h : 


_ . [Dasi-Espuig et al. 

12011 ). The green, blue and red colour 
tones are used for the reconstructions based on the 10 Be, 
sunspot and 14 C data, respectively. The black dotted line marks 
the TSI value at modern solar activity minimum conditions ac¬ 
cording to SORCE/TIM measurements. 


ing on the 10 Be data is the TSI reconstruction bv lshapiro et al 


d201ll) . The magnitude of the secular change in this model 
comes from the difference between semi-empirical model at¬ 
mospheres describing the darkest parts of the solar surface 
and the average quiet Sun (during modern times). The secu¬ 
lar change thus obtained is significantly higher than in all other 
models (3 W m -2 even after the re-assessment by [Judge et al 


2012b . However, the shape of the secular change comes again 
from a linear regression to the 10 Be record, so that it is essen¬ 
tially the same as in the other two models using 10 Be concen¬ 
trations. 

Blue color tones show reconstructions that employ more 
physics-based approaches and are build on the sunspot number. 


Wang et all d2005 ) reconstructed the TSI fro m the magneti c 


flux evolved by a fl ux transport simulation (Sheelev. 2005 ). 
Krivova et al.1 (12010b employed the approach of Solanki et al.1 


( 2002) to compute the magnetic flux from the sun s pot nu m- 
ber and therefrom the technique of iKrivova et al. (12007b to 


compute the TSI and the spectral irradiance. 


(2014) used a flux transport model (Jiang et al 


Dasi-Espu i g et aL 

201 rt EM 


to simulate solar magnetograms to which they applied the 
SATIRE-S model ( Fligge et al. . 2000l:IKrivova et alfl2003l) . 
Finally, Vieira et al. d201 lb adapted the models by 
( 2002b and Krivova et al. ( 2007k for use with 14 C 


Solanki et al 


data (red color tones). The red lines show reconstructions by 
Vieira et al. ( 201 1) based on two different models of the geo¬ 


magnetic field. They hardly diverge over the considered period. 

The magnitude of the secular increase in the TSI since the 
MM differs by roughly 0.5 W m -2 between the models shown 
in the figure. The full range of this variation is actually con¬ 
siderably larger, since the change obtained by [Shapiro et al.1 


2011L in the correction proposed by Judge et al. 2012, is about 
3 W m -2 . More important for the purposes of the present pa¬ 
per is that despite this quantitative difference, the trends shown 
by all reconstructions are qualitatively similar. In other words, 
quite irrespective of the data or the technique used for the 
TSI reconstruction, the TSI is always lower during the MM 
than during the Dalton minimumby on average 0.2-0.3 Wm~ 2 
(more for the model bv lshapiro et al.ll201 1 ). This appears to be 
a rather robust feature. If the sunspot number during the MM 
were as high as proposed by ZP15 then this difference would 
vanish. This contradicts all TSI reconstructions, including ones 
that are not built on sunspot numbers at all, although the un¬ 
certainties in the TSI reconstructions are sufficiently large that 
this cannot be judged to be a very stong constraint (and is not 
entirely independent of the other arguments provided in this 
paper). 


5. Conclusions 

We have revisited the level of solar activity during the Maunder 
minimum, using all the existing, both direct and indirect, 
datasets and evidence to show that the activity was very low, 
significantly lower than during the Dalton minimum or the cur¬ 
rent weak solar cycle # 24. We have confronted the data avail¬ 
able with two scenarios of the solar activity level during the 
MM - the low activity (L-scenario) and the high activity (H- 
scenario, see Sect. [2]). 

The results can be summarized as follows: 

We have evaluated (Sect. 12.IE using a conservative ap¬ 
proach, the fraction of sunspot active days during the MM. The 
fraction appeared small, implying a very low level of sunspot 
activity. 

We have revisited (Sect. 12.21) the telescopic solar observa¬ 
tions during the MM and we conclude that the astronomers of 
the 17th century, especially in its second half, were very un¬ 
likely to be influenced by the religious or philosophical dog¬ 
mas. This is contrary to the claim by ZP15. 

We have discussed that the short gaps in the HS98 database, 
which were interpreted by ZP15 as deliberate omission of 
sunspot records by the 17th century observers for non-scientific 
reasons, are caused by a technical artefact of the database com¬ 
pilation and do not correspond to observational lacunas. 

We have pointed out (Sect. El outdated and erroneous 
information and serious methodological flaws in the analysis 
done by ZP15 that led them to severely overestimate the solar 
activity level during the MM. 

The latitudinal extent of sunspot formation (Sect. E3 also 
points to a very low solar activity during the MM. 

We have shown (Sect. 12.61 ) that East-Asian naked-eye 
sunspot observations cannot be used to assess the exact level 
of solar activity. The existing data do not contradict to the very 
low activity level during the MM. 

We have presented a re-analysis of several documented 
sources, some new to science so far, on the occurrence of lower- 
latitude aurorae on Earth. We have demonstrated that the MM 
indeed displayed very low activity also in terms of auroral 
sightings (Sect. 13.11) compared to normal periods and also to 
the Dalton minimum. 
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We have compared (Sect. 13.41) the estimated heliospheric 
conditions (Sect. 13.31) for the two scenarios, with the actually 
measured cosmogenic isotope data for the period around the 
MM. The comparison is fully consistent with the L-scenario 
but rejects the H-scenario at a very high confidence level. 

We have argued (Sect. 14.11) that the observational facts (very 
low sunspot activity, hemispheric asymmetry of sunspot for¬ 
mation, unusual differential rotation of the solar surface and 
the lack of the K-corona) imply a special mode of the solar 
dynamo during the MM, and disfavor an interpretation of the 
latter as a regular minimum of the centennial Gleissberg cycle 
as proposed by ZP15. 

We have discussed (Sect. 14.21) consequences of the MM for 
the solar irradiance variability which is a crucial point for the 
assessment of s olar variability influence on both globa l and re¬ 
gional climate (lLockwoodll2012l:ISolanki et all. 1 201 3l) . 

We conclude, after careful revision of all the presently 
available datasets for the Maunder minimum, that solar activ¬ 
ity was indeed at an exceptionally low level during that pe¬ 
riod, corresponding to a special Grand minimum mode of solar 
dynamo. The suggestions of a moderate-to-high level of solar 
activity during the Maunder minimum are rejected at a high 
significance level. 
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